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ABSTRACT 

The effect of injection of up to ten percent by volume air in damping 
the amplitude of water borne pressure pulses was studied. Pressure pulses 
were generated using a water hammer model. The water hammer model con- 
sisted of a length of pipe with water flowing through it and a quick clos- 
ing valve located downstream. Measured amounts of air were injected into 
a water stream and the valve was quickly closed. The pressure-time history 
was recorded on oscilloscope photographs at two locations along the pipe. 
A mathematical model was developed for this two phase model. Small amounts 
of air greatly reduced the amplitude of the pressure pulses. By injection 
of one percent by volume of air into the model at 15.5 psia, the pressure 
amplitude was reduced to 8.6 percent cf the theoretical no air amplitude. 


The effect decreased with increased pressure of the system. 
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Lt. Introduction. 

In recent years the United States Navy has been doing a great deal 
of research and development in methods of reducing the noise level of 
machinery for submarines and for anti-submarine ships. The fluid systems 
of these ships transmit noises via pressure pulses in the systems’ con- 
tents. A similar problem concerning fluid borne pressure pulses is that 
of vibrations set up by water hammer in piping systems. The objective of 
this study was to quantitatively determine the effect of air in damping 


water borne pressure pulses. 
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2% Theory. 

To study the effect of air in damping water borne pressure pulses, 
ist: ee to generate a pressure pulse in the liquid. For this 
study, the water hammer model was used to generate the pressure pulse. 
The water hammer model shown in Fig. 1 consists essentially of a head 
tank, a length of pipe, and a quick closing valve. The water is allowed 
to flow through the pipe with a velocity head. If the valve is quickly 
closed, this velocity head is immediately changed to pressure head which 
in turn sends a pressure pulse down the pipe. This pressure wave travels 
down the pipe at the speed of sound in the fluid contained in the system. 

Considering the system shown in Fig. 1, the maximum pressure rise by 


instantaneously closing the valve is given by the following equation, 


Vo Q 
= 1 
Where: iy = Head rise above the initial static no flow head, ft. of water 
lL. = Velocity of the water in the pipe before valve closure, ft/sec. 
S20 Proportionality constant in Newton's second law of motion, 


(1bm- £t)/(1bf-sec2). 

a = Velocity of the pressure pulse, ft/sec. 
The fer ieeeion os this equation as derived by Stepanoff /13/ is included 
in Appendix B. This equation does not include frictional effects. If 
the valve is closed instantaneously, the maximum pressure given by equa- 
tion (1) will be the same along the whole length of the pipe. Obviously, 
a valve cannot be closed instantaneously. However, this same maximum 
pressure rise at a given location A along the pipe will be attained if 
the valve is closed in less time than it takes the pressure pulse to travel 


from position A to the tank and to reflect back to position A. At the 





valve, the time for the pressure pulse to reflect back is 2L/a. Where 
L is the length of the pipe from the valve to the tank. At the center 
of the pipe, the time for the pressure pulse to reflect back is L/a 
Next to the tank the maximum pressure will not be attained unless the 
valve is closed instantaneously. For the maximum pressure given by 
equation (1) to be attained at a given point located at a distance X 
from the valve, the closure time must be less than 2(L-x)/a ° 


The speed of the pressure wave in a pure water pipe system is given 


by 
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QQ = (2) 


Sige, ee 


Water density, slugs/ft3 


Where: ie 


Water bulk modulus, eS . 


K = 
WwW 

D = Inside pipe diameter, in. 

e = Pipe wall thickness, in. 

E = Pipe modulus of elasticity. 


This equation has been derived in Appendix B. This equation can be ex- 
tended to a two phase air-water system. The speed of a pressure pulse 
based upon an isothermal compression of the air by the pressure pulse is 


given by 


- “v0 
= ie + (I= x + mex [ae fa Kale ge lo-ropatxpa (3) 
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Where: x = Volume fraction air at water inlet conditions. 

a, = Velocity of the pressure pulse in the air-water mixture 
with isothermal compression of the air by the pressure 
pulse, ft/sec. 

ave= Velocity of sound in air at the temperature of the fluid, 


ft/sec. 





a_= Velocity of sound in water at the temperature of the fluid, 
ft/sec. 


k = Ratio of specific heats at constant pressure to that at 
constant volume for air. 


Paps = Absolute pressure of the system, Tee ee 
Pa = Air density, slugs/ft> 


ae = Compressibility of water, £t7/1bE. 


The speed of the pressure pulse based upon an adiabatic compression of the 


air by the pressure pulse is given by: 


a;? = @ + (oO sxe fy oe K.,|+ Bd») pu. +XPa | (3a) 


Equations (3) and (3a) have been derived in Appendix B. It is noted that 
equations (3) and (3a) reduce to equation (2) when no air is present. 
Substituting equation (2) into equation (1) and multiplying by the specific 
weight of water CS @)3 gives the maximum pressure rise upon valve closure 


with no air in the pipe. 


Vo Yous 


ag, STs = Jo Pw (A. + ie) 


The maximum pressure rise for an air-water mixture based upon isothermal 





(4) 


compression of air is given by: 


‘fp 


P= Wied (5) 


co}. 


The term ay is given by equation (3). The derivation of eauation (5) is 


in Appendix B. The same equation based upon adiabatic compression of the 


4 


= so ta 








= ¢ eo —- eae | 


- a 
ee ee oe FP 


oe = oe —_—_——-— =_a«'); ——= am & 
<_< <= >_>) a 


(eee @ a: epmrais 














-— Ti ty wl 
— ~ a * 
ps oe 5° ©. =p a> 
SSS ee @ ON ee ee 
-—=—— - 
ee = 
— — —_— Se => ee, 
—_ 4% 





air is given by: 
p _ Ve %eO 
Go 


is given by equation (3a). The derivation of equation (5a) 


(5a) 


The term ay 


is in Appendix B. The ratio of equation (5) to equation (4) gives: 


ee oa xs oe 
ee, w mE s=) 2 
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The ratio of equation (5a) to equation (4) gives: 





si! CO 
Fw ic a 2} 


(6a) 





If equation (6a) is used, adiabatic conditions are assumed or if equation 
(6) is used, isothermal conditions are assumed. Normally, sound traveling 
in the atmosphere is adiabatic. Karplus /4/ conducted experiments with 
air bubbles approximately 0.01 centimeters in diameter in an air-water 
mixture and found that his experimental data was best correlated by assum- 
ing an isothermal compression. The method used korameee air for this 
investigation did not give the small bubble distribution Karplus obtained. 
In Karplus' method, there was a high heat transfer surface per unit bubble 
volume. Therefore, the compression was nearly isothermal in the water. 
With larger bubbles there is less surface per unit bubble volume and the 
situation is probably best described by the adiabatic equation (6a). 
Equation (6a) provides the fundamental theory for the damping of 
pressure pulses in water by air injecticn. By waeeie ine ‘the air injection 
rate, the water flow rate, and the pressure-time history at a given loca- 
tion in the system of Fig. 1, the effect of air in damping a pressure 
Butie in water is obtained. In this report, the experimental pressures 


obtained are compared with the theoretical relationship of equation (6a). 
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3a Equipment. 

To study the effect of air in damping water borne pressure pulses, 
and to verify the applicability of equation (6a) to the water hammer 
model, the system shown in Fig. 1 was constructed. 

The system consisted of 1) a head tank, a pipe and a weigh tank 
system 2) an air injection system 3) a vaive and closure system and 
4) a pressure recording system. A general description for each part of the 
system will be presented here. Exact detailed equipment specifications are 
included in Appendix C. 

The head tank was a tank two feet in diameter with four constant head 
overflow sections. These overflow sections were located nominally at two, 
four, six and eight feet above the main tank outlet line. The main pipe 
length consisted of 59 feet 6 and 7/8 inches of two inch extra strong 
black steel pipe. This pipe was supported at ten feet intervals by a timber 
structure. At the main pipe discharge, a 20 inch by 28 inch by 20 inch 
weigh tank, platform scale and electric timer were used to measure the 
weight flow rate of the water. Water for the system was supplied from the 
city water supply. 

The air injection system contained a Fischer and Porter Company roto- 
meter with a glass ball float for air flow measurement. This rotometer 
had a maximum capacity of 0.325 cubic feet per minute at 14.7 psia and 
70°F. Compressed air from a nominal 100 psig compressor was reduced by a 
regulating valve to approximately fifteen psig. A mercury manometer was 
used to measure the air pressure, and the air temperature was measured 
with a copper-constantan thermocouple junction upstream of the rotometer. 
The air was transferred from the regulating valve to the point of injection 


through one quarter inch vinyl tubing. The air was regulated by a manual 





control valve located immediately following the rotometer. Fig. 1 shows 

a diagram of the system, Fig. 2 shows a photograph of the component parts 
of the air injection system. The air was injected into the pipe near the 
tank through a system of three circles of small holes drilled equally 
spaced around the circumference of the pipe. Details are included in 
Appendix C. A concentric plexiglass cylinder with "0" ring seals surrounds 
these holes. The system was designed so that any one, any two, or any three 
of the circles of holes could be used to inject air at any one time. Fig. 
5 shows a cross sectional view of the air inject annular ring. This method 
for air injection is very simple and could easily be adapted in practical 
applications. 

If the same maximum pressure (neglecting friction losses) was to be 
obtained along one half the length of the pipe, it was necessary to close 
the valve in less than one half the pipe period (<=). This was obtained 
by the use of a modified two-inch Crane quick-closing valve. A steel cy- 
linder two inches © diameter and one inch high was screwed onto the top 
of the stem, The valve was closed by dropping a seven pound weight from 
nine feet onto the cylinder. A spring loaded catch was provided to keep 
the valve closed when the valve was in the shut position. An isolated 
stop stand was provided for stopping the valve. This stop stand minimized 
the effect of the valve stopping impact from being transmitted to the pipe 
and to the valve body. A quick release mechanism was provided to release 
the weight which was guided down through a perforated aluminum tube, one 
and one-quarter inches in diameter. Fig. 3 shows a photograph of the 
valve and closure system. In Appendix G, calculations are provided which 
were used to estimate the closure cine of the valve. The closure time of 
the valve was measured with an oscilloscope and a Polaroid camera. The 
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time sweep was triggered with a photoelectric cell. The closure was 

marked by a 45 volt direct current signal on the vertical plates of the 
oscilloscope. The distance that the sweep moved on the oscilloscope 

from trigger to closure, was photographed with 3000 speed Polaroid film. 

The closure time for the seven pound weight from nine feet was 0.007 + 0.003 
seconds. This agrees closely with the values calculated in Appendix G. 

A Kistler model 401 quartz pressure transducer was installed at three 
locations along the pipe. A Kistler 407-5 low noise cable was used to 
connect the transducer to the Kistler model 651 piezo-calibrator amplifier. 
A shielded cable was used to connect the output of the piezo-calibrator 
amplifier to a Hewlett Packard model 130 oscilloscope. An oscilloscope 
Polaroid camera was used with Polaroid 3000 speed film to take a time 
exposure of the pressure-time history for each run. When the valve was 
completely closed, a 45 volt direct current circuit was completed. This 
triggered the oscilloscope time sweep. Fig. 4 is a photograph of the com- 


ponent parts of the pressure recording system. 
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4, Procedure. 

The objective of the experimental part of this work was to verify 
the applicability of the theoretical model. 

The pressure-time history was recorded at three locations along the 
pipe. Fig. 13 shows these locations. They were numbered starting from 
the valve, one, two, and three respectively. The pressure-time history 
was recorded at each transducer location for five different air-water 
mixtures at each nominal head in the head tank of two, four, six and 
eight feet. Since the effect of small amounts of air was of interest, 
air rates were chosen such that the volume percent air in the system was 
always less than ten percent. 

Because of the effect of as little as 0.1 percent by volume air in 
the water for the no air runs, it was difficult to get good results with 
water taken directly from the water main. For these runs, it was found 
that better results were obtained if the water was left in the tank over- 
night. 

The transducer was calibrated using an Ashcroft dead-weight tester. 
The calibration curves are given in Appendix C, Fig. 15 and 16. 

A calibration curve for the rotometer was available from the manu- 
facturer. Four points on this curve were eneerea with a positive displace- 
ment wet gas meter. The points checked closely, and the curve was used. 

The detailed procedure for each run is given in Appendix D and a 
sample data sheet is shown in Appendix E. 

Only one pressure transducer was used. This transducer was moved 


to the desired location for each run, 
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5 Results. 

The results have been presented in the form of tables, photographs, 
and curves. In Table 1 "transducer locations’ are numbered as shown in 
Appendix A, Fig. 13. The "head" refers to the head of water in the head 
tank. The values given are nominal, and the exact values of the heads are 
shown in Appendix A, Fig. 13. The term Y} is the volume percent of air 
at the water inlet conditions of temperature and pressure. The term Mop 
is the weight percent of air in the mixture. P is the maximum pressure 
rise above the static pressure with no flow in the pipe. The ratio P/P oe 
is the ratio of the maximum pressure obtained when the valve is closed to 
the theoretical maximum pressure with no air, equation (4). 

Results are tabulated for transducer locations one and two. Meaningful 
results were unattainable from location three. The reasons for this will be 
discussed later. 

Fig. 6 shows a series of photographs for nominal heads of two and four 
feet of water and Fig. 7 shows the series of photographs fer nominal heads 
of six and eight feet of water. Each row of photographs is for a given 
head. In the first photograph no air was injected into the water. The 
amount of air increases moving from left to right. The camera used to take 
the photographs takes the pictures such that the sweep of the oscilloscope 
moves from right to left rather than left to right as the oscilloscope 
sweep actually moves. The initial point in Fig. 6 and 7 has been marked by 
lines. The horizontal arrows point in the direction that the sweep moves 
along the horizontal scale. The vertical arrows point in the direction of 
increasing pressure. In the photographs, it should be noted that the pres- 
sure and time scales are different when air is injected from those when no 
air is injected. | 


10 















——_——_ «= = qa -—aae 
ms — > ae aa™ - 
mee ss “Se — @ ae 
¢ & ——— a ae ; 
_— S$. << ax» >i! 
_— > - —_ eae = ee 
ee eee FS | aes Fe Se 
— Ds ht a <a: az 
ck <= Gee Gea aa = 
—— - a 

—~ =i > = ja 


o=_r—n_-<—— = an” 








In Fig. 8 through 11, the volume percent of air was based upon 
the water conditions of temperature and pressure at the inlet to the main 
run of pipe. Fig. 8, 9, 10, and 11 are respectively for nominal heads of 
two, four, six and eight feet of water. Both the theoretical adiabatic 
(Equation (6a)) and theoretical isothermal (Equation (6)) curves have been 
drawn in. The experimental points for transducer locations one and two are 


plotted and marked with symbols. 
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Gr Discussion of Results. 

Transducer location three could not be used to obtain meaningful 
data because the valve could not close fast enough. The valve would 
have had to close in less than 0.0015 seconds to get the maximum pressure 
without interference of the reflected pulse returning from the tank. Hence, 
only data from transducer locations one and two are presented. This sub- 
ject was previously discussed in the theory section. 

As the tabulated results show, two runs were made for each set of run 
conditions. This was done as a means of checking the accuracy and consist- 
ency of the results. 

The weight percent of air in the two phase water-air mixture is includ- 
ed in the tabulated results to show the small weight fraction of air pre- 
sent in the air-water mixtures. 

It was impossible to obtain meaningful results at transducer location 
one with no air in the water. The high frequency pressure pulses set up 
by the valve stem and gate vibrating made it impossible to distinguish 
between the pulse created by stopping the water and that created by the 
vibrating valve parts. With a small amount of air, these high frequency 
pulses were damped out and meaningful results were cbtained. The effect 
of air in damping the high frequency noise is seen in Fig. 6 and 7. 

The pressure-time photographs of Fig.6 and 7 show how the pressure- 
time history at transducer location two appeared on the oscilloscope. The 
photographs show that in the region observed; 1) the greater the percent 
of air present in the mixture, the smaller the pressure amplitude and the 
slower the pressure pulse moves down the pipe. The slowing rate of the 
pressure pulse is shown by the increasing time duration of the pulse as 
the air percent increases. 2) Small amounts of air have a large effect in 
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reducing the amplitude of the pressure pulse and in decreasing the pressure 
pulse speed. 

The photographs taken at transducer location one show the same effects 
as the photographs taken at transducer location two. As a result, the 
photographs taken at transducer location one are not included in this re- 
port. Fig. 6 and 7 show the photographs taken at transducer location two. 
The results taken from the photographs are tabulated in Table l. 

In Fig. 8 through 11, the theoretical curves include no frictional 
losses. Because there are no frictional losses included in the theoretical 
curves, it would be expected that the experimental results should fall 
slightly below the theoretical curve. This is because the frictional 
losses for a two phase mixture are greater t*.an those for a single phase. 
Martinelli experimentally showed this. /10/ As discussed in the theory 
section, the compression of the air by the pressure pulse was probably best 
described as adiabatic because the method of injection allowed for large 
bubbles which gave a small heat transfer area per unit volume. The pres- 
sure pulse travels so fast that the air is compressed for a very short 
time, consequently this permits very little time for heat transfer. 

For the runs when no air was injected into the water, values of 81.7, 
71.7, 74.7 and 79.9 percent of the theoretical maximum pressure were obtain- 
ed at transducer location two for nominal heads of two, four, six, and 
eight feet respectively. Entrained air in the water of about 0.02 percent 
by volume would cause this effect. The entrained air was under pressure 
of about 100 psia in the United States Naval Postgraduate School system. 
When the water pressure was reduced, the air expanded and the volume per- 
cent increased. A glass of water will have air bubbles on the glass sur- 
face if allowed to sit overnight. 
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By allowing the water to stand about fifteen hours in the tank, the 
maximum pressure obtained when no air was injected increased from approxi- 
mately 50 percent to about 80 percent of the theoretical maximum. This 
indicates that the runs in which no air was injected had enough entrained 
air in the water that the exact air content was uncertain. A factor of 
0.02 percent by volume makes a large difference in the maximum pressure 
obtained in the region approaching zero percent of air. This same per- 
cent of uncertainty makes very little difference in the maximum pressure 
in the region greater than about 0.5 percent of air. Therefore, this 
entrained air did not effect the results appreciably for values greater 
than 0.5 percent of air. 

The curves of Fig. 8 through 12 were plotted on a scale in which the 
low ranges of pressure and volume percent were expanded. This was done to 
better show the region of greatest change. A log scale was used in which 
for the ordinate In ( ua + 0.1) was plotted and for the abscissa 
In ory 1.0) was plotted. This was done so that it would be possible to 
have zeros for P/P,_, and nN cn the scales. The results indicate that the 


NT 


points fall close to the theoretical curves in the regions greater than two 
percent of air. In the region of less than two percent of air, the ee 
obtained are more erratic. The uncertainty in the region of 0.5 percent 

by volume air was plus or minus about ten percent in the volume percent 

term, and about plus or minus 11.0 percent in the ee term, Even consider- 
ing these uncertainties, several of the points in the region less than two 
percent by volume air do not fall on the theoretical curve. In the region 

of greater than two percent by volume air the consistency of the data vas 
much better. At four percent by volume air, the uncertainty of the volume 


percent term (Y] ) was about plus or minus 0.8 percent and that of the P/ Pie 
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term was about plus or minus 17.0 percent. These uncertainties are 
based upon the uncertainty analysis of runs 3-9-2 and 3-9-10 given in 
Appendix H. 

Fig. 12 shows the effect of air in damping pressure pulses in systems 
at pressure up to 500 psia. This figure shows curves of (P/P ne V8 Nl 
for pressures of 500, 400, 300, 200, 100, ard 15.5 psia. The volume per- 
cent of air at inlet conditions varies from zero to ten percent. These 
curves show that the lower the pressure in which the pressure pulse is 
being transmitted, the greater is the effect of air in reducing the pressure 
amplitude. As an example, at one percent by volume of air at 15.5 psia, 
the pressure would be reduced to 8.6 percent of the original pressure, 

At 100, 200, 300, 400, and 500 psi the pressure amplitude reduction would 
be to 20.8, 28.7, 34.4, 38.9, and 42.5 percent respectively of the original 
pressure amplitude. 

As an example of the usefulness of air in reducing pressure pulses 
and also to get an idea of the amount of air that would be required, con- 
sider the following example. 

Example 1. Pure water is flowing through a two inch extra strong 

pipe at the rate of 500 gpm. The system has an operating pressure 

of 100 psia and an operating temperature cf 70°F. Fluid borne noises 

are being transmitted in this system. [It is desired to reduce the 

amplitude of these pressure pulses to 20 percent of the original 
amplitude, 

Fig. 12 shows that for the pressure amplitude to be 20 percent 
of the original amplitude, a 1.05 percent by volume mixture must be 


made by the injection of air. 
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00 
+48 qAl/tE = 66.84 ft/min water 


Let x equal the air required in fe Fein at 100 psia and 70°F. 


Then 


100X 
X + 66.84 


~~ — 0.7093 fe i at 100 psia and 70°F 


0.0105 = 


This equals 4.825 feat at 14.7 psia and 70°F. 


The weight of air required would be 0.3628 1bm/min. 
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7. Conclusions. 
Based upon the experimental work and development of equations (6) 
and (6a) the following conclusions have been made. 

1. Small amounts of air in a water system greatly reduce the 
amplitude of a pressure pulse transmitted in the water. Ina 15.5 psia 
water system with one percent by volume air, the pressure is reduced to 
8.6 percent of the theoretical no-air amplitude. 

2. The effect of air in damping pressure pulses decreases as 
the pressure of the system increases as demonstrated by Fig. 12. 

3. There was a reasonable correlation of experimental results 


with theory as shown in Fig. 8 through 11. 
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KEY TO SYSTEM DIAGRAM 
Water storage tank 
Two inch water fill line 
Storage tank overflow piping system 
Overflow drain tank 
Air injection system 
Two inch pipe coupling 
Fischer and Porter rotometer 
Potentiometer reference ice bath 
Air manifold 
Inlet air mercury manometer 
Thermocouple potentiometer 
Kistler Model 401 piezoelectric pressure transducer 
Kistler #470 shielded cable 
Kistler Model #651 piezo-calibrator 
Shielded calibrator-oscilloscope cable 
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Dayton scales 
Weigh tank 
Air inlet themocouple 
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Water inlet Weston thermometer 

Standard ASA two inch extra strong steel pipe 
One quarter inch vinyl tubing 

Two inch gate valves 

Air regulating valve 

Quick closing air stop valve 

Air flow regulating valve 

Water inlet regulating valve 

Weigh tank two inch discharge gate valves 


Oscilloscope sweep trigger 
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KEY TO EQUIPMENT PHOTOGRAPHS 
Description 
Thermocouple potentiometer 
Fisher and Porter rotometer 
Potentiometer reference ice bath 
Inlet air mercury manometer 
Air injection system 
Hewlett Packard model 130 oscilloscope 
Kistler model #651 piezo-calibrator 
Kistler #470 shielded cable 
Kistler model 401 piezo-electric pressure transducer 
Modified Crane quick closing valve 
Falling weight gvide tube with weight in it 
Weigh tank 
Dayton scales 
Weigh tank discharge valve 


Electric timer 
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Table 1. 


Run No. 


3-21-32 
3-21-33 
3-21-34 
3-21-35 
3-21-36 
3-21-37 
3-21-38 
3-21-39 
3-21-40 
3-21-41 
3-21-18 
3-21-19 
3-21-20 
3-21-21 
3-21-22 
3-21-23 
3-21-25 
3-21-26 
3-21-27 
3=21-3 

3-21-4 

3-21-5 


3-21-6 


Tabulated Results 


Trans- 
ducer 
Location 
1 


1 


Head 
tte 


Yt 


1309 


1.09 


5.14 


7219 
9.23 
9728 
0.664 
0.664 
2.02 


2501 


4.77 


6.38 


0.538 
0.536 
1,53 


io8 


a3 


10.05 
10.04 
13.16 
13226 
183 


T363 


6.85 
9.3) 
2e31 
0.774 
0.769 
228 


2.28 


10.6 
10.6 


10.6 


S056 
28.8 
26.4 


24.6 


P/P 


0.0498 
0.0498 
0.0498 
0.0526 
0.0413 
0.0385 
0.0343 
0.0343 
0.0300 
0.0300 
0.0704 
0.0704 
0.0462 
0.0497 
0.0349 
0.0374 
0.0327 | 
07,0327. 
0.0327 
0.0765 
0.0720 
0.0660 


0.0615 


rable 1. 
Run No. Trans- 
ducer 
Location 
3-21-7 1 
3-21-8 1 
3-21-13 1 
3-21-14 1 
3-21-15 1 
3-21-16 yi 
3-13-3 i 
3-13-4 i 
3-13-5 1 
3-13-6 1 
3-13-7 i 
3-13-8 1 
Bo-9 1 
3-13-10 1 
3-13-11 1 
3-13-12 1 
2-27-4 2 
4-1-4 2 
4-1-5 2 
2-28-4 z 
2-28-5 2 
2-28-13 2 
2-28-14 2 


Tabulated Results (continued) 


% 


2,00 
2.60 
3.68 
3.68 


4.76 


0.475 


0.475 


28 


2.14 
S70!) 
3.01 


e791 


0,294 
20 
el 


oro 


34 


in 


%x10° 3 


3.84 
3.84 
5.49 


5.49 


8.07 
O77 31 
0.731 


2.04 


300 
330 


4.77 


Ged 


On 352 
0.382 
E65 
1.58 


4,25 


18.0 
15.6 
1 a0 
1250 
14.4 
55.6 
41.7 


So0 


25d 
Zo a3 


roe 


eed 
de 
174.0 
3.41. 


45.4 


Ta78 
S207 


997 


Pye 


0.0450 
0.0450 
0.0390 
0.0390 
0.0300 
0.0360 
0.0773 
0.0901 
0.0724 
0.0747 
0.0555 
0.0633 
0.0413 
0.0451 
0.0387 
0.0413 
07617 

0.202 

0.213 

0.0410 
0.0365 
0.0379 


0.0468 





Table 1 Tabulated Results (continued) 


Run No. 


2-28-15 
2-28-16 
2-28-17 
2-28-18 
2-28-19 
3-27-3 


4-1-6 


2-28-26 
2-28-28 
2-28-29 
2-7 6=30 
2-28-31 
2-28-32 
2-28-33 
2-28-34 
2-28-25 


3-27-2 


3-3-6 
3-3-7 


3-3-8 


Trans- 
ducer 
Location 


2 


2 


Head 
fe. 


YN 
%, 


520 
5.18 
7.20 
Teo 


9.54 


0.277 
One) 
0.769 
27087 
2,067 
3.49 

Bro 


5SmOZ 


6.24 


6.20 


0.136 
0.596 


0.602 


WE 
*%X1LU 


Teco 
Tee 
LOL 2 
TOF1 


ae 


O57 


Orza7 


5320 
TeoD 
7.48 
9.21 


907 


07207 
0.877 


Ofo73 


2.41 


3.90 


228 

LOS. 
84.1 
37 39 
1969 
13.7 
Te 


12.7 


rOe9 
2281 
161.8 
28.6 
50.6 
16.9 
20,52 


22.0 


Ee 


0.0266 
0.0294 
O70283 
OFOZ52 


0.0280 


Ost 

0.261 

0.118 

0.0626 
0.0494 
0.0380 
0.0399 
0.0286 
0.0380 
0.0286 
0.0343 
OR TAR) 

0.404 

0.0715 
0.0956 
0.0422 
0.0655 


0.0550 


Table 1 Tabulated Results (continued) 


Run No. nae ie An Vy ws P ae 
Location 
3-3-9 2 6 2.96 3.89 16.0 0.0400 
3-3-10 2 6 3.61 ao 17.8 0.0445 
3-3-11 2 6 3.61 5.91 126 0.0445 
3-3-12 2 6 4.70 7.24 15.4 0.0385 
3-3-13 2 6 4.76 Tae 10.6 0.0265 
3-27-1 2 8 0 0 370 0.799 
4-1-10 a 8 O. 117 02179 356.9 On77 1 
3-9-2 2 8 0.485 Os750 o149 O.112 
3-9-3 2 8 0.484 0.747 23 0.0913 
3-9-5 2 8 1.26 Zod 48.9 0.106 
3-9-6 2 8 2e13 3.41 26.7 0.0576 
3-9-7 2 8 2.15 3.45 27.9 0.0602 
3-9-8 2 8 3.02 4.76 192 0.0421 
3 9-9 2 8 5.0 4.76 20.3 0.0460 
3-9-10 2 8 3.89 6.37 17 el 0.0462 
3-9-11 2 8 sooo 6.37 ai.) 0.0462 
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APPENDIX B 
THEORY DERIVATIONS 
1. The following is the derivation for equation (1) of the theory 


section as given by Stepanoff /13/. Newton's second law is as follows: 


r= dmv) (B-1) 


ae 


Gite 
Where: F = force, Ibt. 
V = Velocity, ft/sec. 
t, = Time, sec. 


M = Mass, slugs. 


The velocity through the pipe we) before valve closure is constant. The 
force (increase in pressure at the valve) is equal to the time rate change 
of momentum of the moving mass. The pressure-force relationship is as 


follows: 


F=AH, (B-2) 


At time ts after closure the portion of the fluid in the length of pipe, X, 


will have its velocity reduced to zero, then: 


Fenths AEUE) 


Where: A Pipe internal cross sectional area, fee 


ox 
il 


Specific weight of fluid flowing, Tbm/ ft. 
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H = Head rise above initial no flow static head, Ho» ft of 
Y water 


g_ = Proportionality constant in Newton's second law of motion, 
(lbm- ft)/ (1bf£-sec2) 


X = Distance to any position on the pipe measured from the 
valve, ft. 





ty = Time after closure, sec. 
ve = Velocity of fluid in the pipe before valve closure, ft/sec. 
The term a in equation (B-3) is the mass flow rate. In equation 
(B-3), x/t, is the velocity of the pressure pulse, a. Substituting a and 
cancelling the Y¥ A term on both sides of the equation (B-3) gives: 
Wee oe = (1) 


2. The derivation of equation (2) in the theory section as given by 
Stepanoff /13/, is as follows, After instantaneous valve closure a region 
of compressed water is created adjacent to the valve and extends a distance 


X at time ty after closure. The pipe in this region expands under pres- 


sure from Do to Di: From the instant of closure until ty a volume of wa- 


ter equal to AV t passes through the original cross section at X. There 


1 
is no outlet for the water; therefore, the volume must be accommodated by 


expanding the pipe radially and by compressing the water. 


Incremental i x ——>| 
Increase due 


to com pression 


of original mass 
(Axe | 
K 





Incremental increase 


due to expansion of 
pe volume 


TN ah) 24h 





Original mass 
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The increase in volume caused by the compression of the water in the pipe 


is obtained from the definition of bulk modulus as follows: 


ae a gla’ 
Ky = —Cae = Gw AGy (B-4) 


Where: Ky = Water bulk modulus, ieee 


3 
G, = Volume, ft 


4G, = Incremental change in volume, ft? 


P = Incremental change in pressure caysed by the velocity head 
changing to pressure head, lbf/ft’. 


Solving equation (B-4) for AG, gives: 


L5G5 = <= (B-5) 


but the original volume G, is equal to AX. Substituting AX for G, gives: 


yan oe = -~AX P_ (B-6) 
Kw 


The negative sign indicates that AG decreases with increasing pressure. 
The decreased volume is used to accommodate the additional mass that 
enters the section of Pipe X feet long. The volume made available by the 


compression of the original water is: 


AG. = AXP (B-7) 
Kw 


The increased volume caused by the expanding of the pipe is; 


AGp = £ (D*-p3) X (B-8) 


By equating the additional volume of water that must be accommodated 
(AV ty) to the volume that the pipe expands plus the volume that the 


water compresses the following is obtained: 
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A iy ae 2 
AV,t; = Aan a 4 (0-0. X (B-9) 
w 
For the pipe 
Se 
1 
—————— PO, ee 
i 
Where: S = Pipe hoop stress, lbf/ft'. 
e = Pipe wall thickness, ft. 
Hooke's law for the pipe expansion is: 
S ‘a p,-O 
= p( Soe (B11) 
Ds 
Where: E = Modulus of elasticity, wey ee 
Combining equations (B-10) and (B-11) gives: 
PD, pat D, a Do (B-12) 
2eE io. 


2 2 : 
In equation (B-9) dD, =P. may be written (D,+D,)(D,-D,)- 


Since D,-D, is very small compared to dD, or Do the term (D,+D.) may be 
written 2D: Then the last tem of equation (B-9) becomes; 
Tt A Dela (OD) — age (0 - 0:) (B- 13) 
4 fe 
Multiplying both sides of equation (B-12) by D, gives: 
Bon 
D,- 9, = ae. (B-14) 
2eE 


Substituting equation (B-14) into equation (B-13) and substituting A for 


ie 
the cross sectional area term, m0, gives 


4 
Xe TT ( p;- D2) — A A Do Pp (B-15) 


4. Se 


Substituting equation (B-15) for the last term in equation (B-9) gives: 
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XPA 4 XD.PA 


V,AtL = K., Ze 
Do : 
ee a ee) (B-16) 


The term X/t, is the wave velocity, a. Substituting and solving for a 
gives: 
ae 
Kw CE 
If equation (1) is multiplied by the specific weight of water ( Oe the 


following is obtained: 





Q Ve ye 
— = Ooo (B-18) 
Hy % =P ac 
Solving (B-18) for v gives: 
Y= xigaoes (B-19) 
% ew 
Sw 
Substituting (B-19) into (B-17) and letting D equal »% and v equal Go 
gives: 
i 
C= 
(2) 


Z 
Pole ec 
The maximum theoretical pressure rise, -(P..) is given by substituting equa- « 


tion (2) into equation (B-18). 


NT 9. /p (+ n 2) 
DEN hy e 
The experimental results were expressed as P/P__. 


a oe 
£ = Pp fe Bie E) ae 





The term P is the experimental pressure rise above the no flow static pres- 
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3. The following is the derivation of equations (3), (3a), (5), (5a), 
(6), and (6a) of the theory section. These derivations are based upon 
work done by Karplus /4/. His equations were modified to make them apply 
to the water hammer model. The following is a form of equations (2) for 


an air-water mixture. 
ne D 
a = pel X. + = | (B-21) 


Where: Subscript c represents the air-water mixture. 


i= - = Air-water mixture compressibility, ft" /1be. 
c 
The density of the mixture CP? is given by: 


ie = Mat Mw = Mo ++ Mw (B-22) 
Ge GC. ‘ 


Where: G Volume, ft. 


M = Mass, slugs. 


subscripts 


a = Air 
W = Water 


Because mass is the product of density and Volume pe canes: 


jem a tae 29) 


— 


The volume fraction of air is defined as: 
= Se (B-24) 
The volume fraction water is defined as: 
(4-x) = Sw (B-25) 
Ge 
Substituting (B-24) and (B-25) into (B-23) gives: 
(a X Pp. + O-X) Pw (B-26) 


b4 
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Compressibility is defined as follows: 


i — ~(dG. + dGw) 








GedP G.dP 
= dAGw ($+ = AGa (Es) (B-27) 


Substituting equations (B-24) and (B-25) into (B-27) gives: 


Ko = (1-xX)®, + ne aes) 


Substituting equations (B-26) and (B-28) into equation (B-21) gives: 


anal 1x) Pw - X Pal| x X, + (-x) Kw + # CERe?) 


Rearranging equation (B-29) gives: 


a= x2p. Ky +0-%) Py Nw + X(I-a) [Pa Kw + Pw HID rrpa xp] Cay 


For isothermal conditions and assuming the ideal gas law applies: 


en = 4 7 (B-31) 
Pras 9 
Where: Pans = Absolute pressure of system, l1bf/ft™. 
Meg = Isothermal compressibility of air, fe“ /ibe 


The velocity of sound in air is given as: 


-2 
— B~ 32 
a Kaa Pa ad 
Where XnaX Adiabatic compressibility of ain ft-/ ibe. 


For adiabatic conditions and assuming the ideal gas law applies: 
XX a (B-33) 
“ R Fags 
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Where: k = ratio of specific heat of air at constant pressure to that 
at constant volume, 1.403 for air 


By combining equations (B-31) and (B-33), the following is obtained: 


Kra = Kaa (8-34) 


The velocity of sound in water is given by: 


Pw Ra (B-35) 
For isothermal conditions substitute (B-32), (B-33), (B-34), and (B-35) into 


(B-30). This gives: 
a7 xk + 0" 5 xoxo Le + X,, Ee Ble-» +X | 

= Q Ow Fans fa *\w CE Pw Pa 
Using the adiabatic conditions for equations (B-30) gives the following 


results: 


gee 1-x ee | _ | 
2 x tage fe x( x)| | 8,1 Ps Kw | + =e |G *) Pus +x | (3a) 


For isothermal conditions for an air-water mixture equation (B-18) becomes: 





Where: ay is given by equation (3). 
Substituting equation (B-36) into (B-20) gives: 
fe Ve SeQs i o) 
ck mw /Pu (K+ i 


Where: a, is given by equation (3). 


For adiabatic conditions: 
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CE) = le Gy + ce) . 


Where: a, is given by equation (3a). 
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APPENDIX C 
EQUIPMENT 

1. The head tank consisted of a tank which was two feet in diameter at the 
lower end and one foot ten and one-half inches at the upper end. The upper 
part of the tank was light sheet metal and the lower part of the tank was 
plate steel. Two inch overflow lines were provided at approximately two 
foot intervals as is shown in Appendix A. These overflow lines were stan- 
dard two inch galvanized steel pipe. Each overboard line had a two inch 
gate valve in the line. All the overboard lines ran into a common line which 
discharged into a collecting tank open to the atmosphere. This collecting 
tank had a drain which discharged into the main sewage system. A Weston 
thermometer located near the bottom of the tank, was used to measure the 
inlet water temperature, 

The main run of pipe consisted of two inch extra strong black steel 
pipe. The outsidediameter was 2,375 inches, the inside diameter was 1.939 
inches. The wall thickness was 0.218 inches. Three sections of pipe were 
welded together. Care was taken to insure good butting of the individual 

on 
sections. At the transducer locations, transducer fittings one inch in 
diameter and one inch high were welded into the pipe. These fittings had 
standard 14 mm by 1.25 automotive spark plug threads. 

The water weigh system consisted of a Toledo Scales Company platform 
scales with a capacity of 1600 pounds and a steel plate rectangular tank, 
20 inches by 28 inches by 20 inches. A 3 inch by 3 inch section was cut 
out of the top to permit entrance of the discharge pipe. This tank had 
two, 2 inch gate valves for discharge. They discharged into drain lines to 
the sewage system. 
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The following scale information is provided, 


Manufacturer Toledo Scales Company 
Toledo, Ohio 


Model 31-1811FC 
Serial Number 2180-0-225 
Capacity 0-1600 pounds 


An electric timer, which measured to the nearest one tenth da 
second, was used to measure the time to fill the weigh tank with water. 
2. The air injection system is shown schematically as part of Fig. 1. 
A photograph of the air injection system is shown in Fig. 2. The source of 
air was from the 125 psi air system permanently installed in Building 500, 
at the United States Naval Postgraduate School. This system was reduced 
to approximately fifteen psi by a regulating valve. One-quarter inch vinyl 
tubing was used from the regulating valve to a pipe manifold. A copper- 
constantan thermocouple was installed in the pipe manifold to measure the 
temperature of the air. A Rubicon Company catalog number 2732 potentiometer 
was used to measure the thermocouple output. This is item 1 in Fig. 2. 
~To measure the pressure of the air a 50 inch mercury manometer was used. 
A one-quarter inch vinyl tubing line was tapped into the air manifold. This 
line was connected to the mercury manometer, The other side of the manometer 
was open to the atmosphere. Item 4 of Fig. 2 shows this manometer. From 
the manifold a one-quarter inch vinyl line led to the rotometer. The follow- 


ing data is provided. 


Manufacturer Fisher-Porter Company 
Flowrotor tube number 2-F-1/4-20-5/70 

Float glass 

Range (14.7 psi, 70°F) 0-0. 325 ft>/min. 
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The manufacturer's calibration curve, converted from cubic centimeters per 
minute to cubic feet per minute, was used. Fig. 14 is a copy of this 
curve. Four points on this curve were checked using American Meter Com- 
pany wet gas meter Al-20-11472. 

The flow rate was corrected to 14.7 psia and 70°R for each run using 


the following relationship suggested by the manufacturer /2/., 


2.= 0 / FE) E) c 


fe mea at 70°F, 14.7 psia 


oe 

Ea = pressure at which measurement was made, psia 

t. = temperature at which measurement was made, °R. 

QF = flow rate in ee aie taken from the curve of Fig. 14. 


This curve was entered with the rotometer reading. 

Item number 2 in Fig. 2 shows the rotometer. A one-quarter inch vinyl 
tubing connected the rotometer with a manual control valve. The control 
valve was used to control the air flow rate. From the control valve a 
one-quarter inch vinyl tube led to: manifold section in the pipe. A 
stop valve was provided for each of these separate lines. Item 5 of Fig. 

2 shows the equipment used for the injection of the air into the pipe. 
Three separate rings of holes were used to inject the air. One ring con- 
sisted of two holes, each with a diameter of 0.04 inches. The second ring 
consisted of four holes, each with a diameter of 0.07 inches. The third 
ring consisted of eight holes with diameters of 0.07 inches. These holes 
were equally spaced around the circumference of the pipe. A plexiglass 
ring, four inches in diameter, with a circumferential groove for air injec- 
tion was provided for each ring of holes. These rings were sealed with "0" 


rings. Fig. 5 shows a cross sectional view of this system. The system was 
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designed so that any one of the rings could be used. However, it was 
not necessary to do this. All experimental work was done with air enter- 
ing all three rings of holes. 
3. The valve and closure system is shown in Fig. 3. The quick closing 
valve was a modified standard Crane two inch quick closing gate valve. 
The valve was modified by screwing on a steel disc two inches in diameter 
and one inch high on the top of the stem. This disc was the point of im- 
pact for the falling weight which closed the valve. A mechanism was de- 
signed to hold the valve closed when it was in the shut position. A stand 
that was completely isolated from the piping system was used to stop the 
valve. By using this method, the full impact of the stopping valve was 
not transmitted to the piping system. A seven pound cylindrical steel 
bar one inch in diameter by 24 inchs long was used as the weight for clos- 
ing the valve. This weight was dropped through an aluminum pipe, one 
and one-quarter inches in diameter. This aluminum pipe contained numer- 
ous air holes to allow the weight to fall freely. Appendix G contains pre- 
liminary calculations used to estimate the closure time. The weight was 
dropped from nine feet for all runs. The closure time of the valve was 
measured een an oscilloscope. The time sweep of the oscilloscope was 
triggered by a photoelectric cell. A contact completed a 45 volt direct 
current circuit when the valve was closed. This voltage was put across 
the vertical plates of the oscilloscope. A Polaroid time exposure was 
made of the oscilloscope, the closure time was taken from the photograph. 
Polaroid 3000 speed film was used. 

The quick release mechanism consisted of a steel pin which held the 
weight nine feet above the valve. A line was attached to this pin. When 


it was desired to drop the weight, the pin was pulled and the weight fell. 
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A nylon string was attached to the weight. This was used with a 
pulley to pull the weight to the desired height above the valve. 

4. Selected characteristics of the pressure sensing equipment are given 
below. For complete details refer to manufacturer's instruction book /5/. 


Kistler 401 Pressure Transducer 


pressure range 


0.0-3000 psi 


overload capacity 507 
linearity 1%, 
natural frequency 48,000 cps 
maximum continuous 

temperature 500°F 
maximum intermittent 

temperature 3000 °F 
length 2s06Mitr, 
diameter 0625 sin, 
weight lw oz. 
serial number #1044 


Kistler Model 651 Piezo Calibrator 


input voltage, maximum 0.6 volts 
output voltage, maximum 0.3 volts 
voltage gain OD 

input impedance to anae 
output impedance 107 ohms 
frequency range 0-10 ke 
linearity 1% 
dimensions 7 xol0ex 67 1n, 
weight 4.7 lbm 





Hewlett Packard Model 130 Oscilloscope 
internal sweep accuracy 5% 


“Vertical and Horizontal Amplifiers 


frequency response 0-300 ke 
sensitivity Imv/cm-20 volts/cm 
accuracy 57% 

input impedance 1 megohm 


Cathode Ray Tube 


power requirements 115/230 volts, AC 
50/400 cps, approximately 
175 watts 
dimensions 9° 3/40in, by 15\in. by. 21-1/4 in. 
weight 39 1bm 


Kistler transducer model 401 serial number 1044 was calibrated using 
an Ashcroft dead-weight tester. The system used for calibration was the 
same as shown in Fig. 4 except the transducer was inserted into the 14 mm 
dead-weight adapter fitting instead of the pipe. Sweeney /14/ describes 
the procedure for using the dead-weight tester. The important thing in 
calibrating was to rotate the weights. when loading the transducer. The 
transducer was rapidly loaded so that there would be no error due to 
charge leakage from the crystal. The transducer was calibrated for the 
calibrator range scales of 10 and 10(.fhe 10 range covers 0-125 psig and 
100 range covers 0-1250 psig. Fig. 15 and 16 show these calibration 
curves. These curves show the readings on the calibrator for a given 
pressure. It was necessary to determine the number of calitration divisions 
that correspond to a centimeter deflection on the oscilloscope. Knowing 
this and using the calibration curves, the psig/centimeter on the oscillo- 


scope was calculated. For the oscilloscope used for this experimental 
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work, the following values were obtained and used. These 


Kistler model 401 serial number 1044 pressure transducer. 


Calibrated 
range 


10 


10 


10 


100 


100 


100 


Oscillescope 


scale 


5 mv/em 
10 mv/cm 
20 mv/em 

5 my/em 

10 mv/cm 


20 mv/em 


37 


Calibration 
divisions/ 

10 em oscillo- 
scope deflection 


ae, 
Soet 
W]e 
ye 
34.4 


67.8 


values are for 


psig/cm 


54.4 


106.9 
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APPENDIX D 


DETAILED PROCEDURE FOR A RUN 


The detailed step by step procedure for each run was as follows: 


Mee 


TO; 
Ee. 
125 


13. 


14. 
i 
16. 


ye 


Record atmospheric temperature and pressure. 

Fill the tank to the level desired. 

Open quick closing stop valve. 

Adjust water inlet to the tank until the overflow at the 
desired level is just a trickle. 


Adjust the air rotometer to desired flow rate with the 


control valve immediately following the rotometer in the air 


inlet system. 

Insure weigh tank discharge valves are completely open and 
that the tank is as low as possible. 

Set electric timer to zero. 

Close weigh tank discharge valves. 

Start timer and record initial weight on scales when the 
timer is started. 

Fill weigh tank until it is nearly full. 

Record weight and time at nearly full position. 

Open weigh tank discharge valves. 

Lift weight to nine foot position and insert quick release 
pin. 

Record inlet and outlet water temperatures. 

Record rotometer setting (read ball top). 

Record air mercury manometer readings. 


Record air temperature. 
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Le. 


rg, 


20. 


21: 


Ze. 


23% 


24, 


25, 


267. 
27 


28. 


29. 
50; 


ale 


Set calibrator range to desired position. Insure 
the meter is adjusted to 1.0. Allow a ten minute warm up 
period for the calibrator and the oscilloscope. 

Set oscilloscope on desired vertical scale. 

Adjust oscilloscope sweep rate to the desired time scale. 

Insure the SYNC switch is set to external trigger (EXT) on 
the oscilloscope. 

Insure that battery is properly connected to trigger 
circuit of the oscilloscope, 

Set the horizontal sensitivity of the oscilloscope to 
internal sweep. 

Check system by checking the calibration reading for ten 
centimeters deflection on the vertical scale of the oscillo- 
scope. 

Adjust vertical and horizontal position until the oscillo- 
scope trace is in desired position. 

Set camera on time exposure, 

Open camera shutter. . 

Close valve by dropping the weight with the quick release 
mechanism. 

Close camera shutter. 

Open quick closing valve. 


Develop Polaroid photograph. 


oo 





. . R 
Py |. Pine Sy TELL EL " 
TTT LL LEER YEE yd! { 
4 ; 
© . o °o 
«o e~ aw 







11 
21. 


18.8) 21.4 


ae 

130 | 130 | 110 
ose 
ree | 

sa 
re 
i] 
10 


[a 
ee 
. | lo 
39.1| 391 39.11 39.11 39. 


a 

re 

Lt 

_ 
70.6 | 70.6 | 70.6 | 70.6 | 70.6 

0.1 | 0.1 | 0.2 | 0.2 | 0.2 | 0,2] 

eee 
a 


13.99 
1262 
i 
a 
300 | 
56 
sa 





39.1} 39.1 


esl 
mare 
Ls 
ea 
Ea 
70.4 
14.41 
14.32 
12.2 
7. os 
i | Sal 
130 
300 
18.4| 18.7 
| 56 [sé | 
| 56 [56 | 56 
2 | 2218 5 
ee 
10 | 10 | 10 





SL! 







120 
jas0_| 260. 
“6 
se | 
ae 
Ee 
0 
sea 
10 
a 
_ 


290" | 
18.1 
_ 
mee 
aaa 
a 
ae 





70.6 | 70.6 | 70.6 | 70.6 


r= 

ae 

Ee 

= 

ee a eS gee | 

| 130 [120 |120 | 

| 298 | 

| 17.5 18.1 


14.50| 14.10-13 


290 









ca 
hel | 
| 3961) 3962 
. 0.1 
10 
[ 
a 





UI UD 
= 
28 
ea 
ws 
Fala 

in Hg 

in Hg 

ei 


6 
|e 
EES AS te] 

| ex 


VERTICAL SCALE 





ROTO} 


HORIZONTAL, EXPERIMENTAL DATA SHEET 


ears 
DIK E 
_- 
—T* 
‘= 





(19-8 *A3Y) SOI SdN ONZI 





APPENDIX F 
SAMPLE CALCULATIONS 
Sample calculations have been included for run number 3-9-10. The 


data for this run has been included on the sample data sheet in Appendix 


Ey 
P ATM = 30226441 nae 
tary = 70. 0°F 
c = -0.129 Barometer temperature correction 


Calculate atmospheric pressure in psia 


(P (in Hg) + c)(0.4913) = Pay psia) 


ATM 
(30.264-0.129) (0.4913) = 14.81 psia 

where 0.4913 is the conversion factor for converting pressure in inches 

of Hg to psia. 


Calculate absolute pressure of air. 


(P)a + Poa) (0.4913) + P = P_(psia) 


ATM 


Pia and Poa are the inches of mercury measured on each leg of the air 


manometer. 
Pia S262 “ine he 
Pos = 12.55 sie He 
P ary = 14,81 psia 
Ee = pressure at which air was measured, psia 


(12.62 + 12.55)(0.4913) + 14.81 = 27.18 psia 


Q, = Q, | (=) 
A147 lig 













—_« et ee _oO 
1 a Ts Fe “ey 


/_ Sa... 
~~, = 
& <= Os 
ange 


ie 
Ee a 


where: 
Q. = air flow rate at 14.7 psia and 70°F, ft>/min. 
Q = uncorrected air flow rate taken from Fig. 14, 
Enter Fig. 14 with rotometer setting. 


T, = Absolute temperature at which air flow rate was measured, 








°R 
2 a = oe 442 ft3 
Q, = On325 14.7 "5 30.b =e ft?/min 
Calculate mass flow rate of air 
Ma a Q4 Ms 


where: 


Mass flow rate of air, 1lbm/min 


we 


Specific gravity of air at 14.7 psia and 70°F. 


¥, 
Ma 


(0.442) (0.0752) = 0.03403 1bm/min 


Calculate volume flow rate at inlet water conditions 





4.7 Tw 
Saw = 8a( “BF ) G35) 
oY = w/ \530 
where: 
Oe = Air volume flow rate at conditions of water inlet 
PY = Absolute pressure of water at air inlet, psia 
caer Absolute temperature of air at water inlet, °R 
AAT (oe 3 
Q aw = 0.442 1B 265) S30 = 0.3463 ft” /min 
Calculate the water mass flow rate 
We.. Mw (60 
Ww ° 
° Bi 
; m,, = Water mass flow rate, lbm/min 
",, = Weight of water weighed in time tis lbm 
ty = time to measure M pounds of water, sec 


62 










el) oe 







i‘) a4 aT cae 
—— a Th oie AY 
es 


— <a 2 
e a ee 








— . 
~ 
é 
i ¢ i | al 


= —— mini 


m= (92) (o) = 534 1bm/min 


W Set 


Calculate water volume flow rate 


where: 


Water volume flow rate, fein 


Qy = 

¥ we Specific weight of water, em ee 
534 37a. 

a = C7 = 6.56 ft /min 


Calculate the volume percent air in mixture 


Q.4 6100) 


a rs Qaw Qy 


where: 


al = Volume percent air at water inlet conditions 7% 


eco ae 
Vln= “eSaea 45.56) = 289% 


Calculate the weight percent air in the mixture 


ag - 


Mo,(100) 


We, + mw 


where: 


ne 
r 


Calculate the maximum pressure 


Weight percent air in air water mixture, % 


0.03403) (100 7 
WI (0.03403 + 534) — 0.00637% 


P = deC+-J 

P = Pressure above the static no flow pressure, psi 

C = Psi/cm of oscilloscope deflection, psi/cm 

d = Deflection of oscilloscope, cm 

J = Correction for velocity head at the location when the 


oscilloscope scale was zeroed, psi. 
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To determine the J correction, the difference in static pressure with flow 

and with no flow was experimentally measured for each head at each loca- 

tion. It was necessary to make this correction because the oscilloscope 

was zeroed for each run when there was velocity head in the pipe. 
Peel o.0)-2,1 ="916,5 psi 


Calculate pressure ratio P/P ion 


p p 


P/P = —_—-> = oe 
= Hy bw (Yew ( 4 ) 
WA: ae (Ree) 
Pipe, = ee 
Va. Sty 


where: 
Kb = Velocity of the water through the pipe with no air, 
determined experimentally, ft/sec 
ie = Density of the water, slugs/ft~ 


K,= Bulk modulus of water, 1b/ ee" 


D = Inside diameter of pipe, in. 
e = Pipe wall thickness, in, 
E = Modulus of elasticity of pipe Ib£/£t? 


Hy = Theoretical maximum head rise above the initial static 


head, equation 1, ft of water 


] 1.239 
<)>) ae 1o.5) (32.2 (312,000)(144) _* ©.219)(22.4% 109144) 
es (7. 52) (62.4) 
P/Pn = 0.0445 
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APPENDIX G 
VALVE CLOSURE 
The water hammer model of Fig. 1 was used to generate the pressure 


pulses, For the maximum pressure to be theoretically possible, the model 


required a valve which would close in less than 26L=x) seconds, 
where: L = The length of the pipe from the tank to the valve, ft. 


X 


The distance from a given location to the valve, ft. 
a = The velocity of the pressure pulse in the pipe, ft/sec. 

The speed of the pressure pulse in pure water in a two inch extra strong 
steel pipe is 4580 ft/sec. For the maximum pressure as given by equation 
(1) to be possible at transducer location two, the valve had to close to 
less than 0.0128 seconds. If there was 0.05 percent by volume air at the 
water inlet conditions tte closure time required at transducer location 
two for the maximum pressure to be possible was 0.0351 seconds. With ten 
percent by volume air, the closure time required at transducer location 
two was 0.438 seconds or less. Since the pressure wave travels fastest 
in water with no air, this was the limiting condition. Therefore, to use 
transducer location two, a closure time of less than 0.0128 seconds was 
required, 

A free falling seven pound weight starting with zero initial velocity 


and falling nine feet would have a velocity calculated as follows: 


oe 2gs_ 


Veo= We(32.2)\(2)) = 24,0 8bt, sec 


where g = acceleration of gravity, ft/sec 


s = distance the weight falls, ft 
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V = velocity of the falling weight, ft/sec 
With a 24.0 ft/sec. velocity, the time of closure would be: 


ee ee 
‘ V 


where t. = time for valve gate to travel S$, distance, sec 


Sy 


distance valve gate travels for closure, ft 


V = velocity of valve gate 


i 


t= (1/6) (1/24) = 0.0069 sec, 
Assuming an plastic collision, the following calculations are made. 
MV = (M, + MV, 
M 
1 
Ve 
2 M)+M, il 


Assuming M, = 1.75 1bm = 1. 75/32.2 slugs 


My 7 Wom =7/32.2 Slugs 


3 a) = 19.2 ft/sec. 


(74+1.75) 
The closure time for the valve would be: 


t= S,/V, = = 0.0087 sec. 


i was.2) 
The packing on the valve was loose enough that the valve closed freely 
with very little resistance. The weight of the valve stem and gate was 
1.75 tbm: 

Experimentally the closure time was measured as 0.007 seconds with an 
uncertainty of 0.003 seconds. The method of measuring the closure time 


? 


has been discussed in Appendix C. 
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APPENDIX H 
UNCERTAINTY ANALYSIS 
The uncertainty for runs 3-9-2 and 3-9-10 are included here as re- 
presentative examples of the uncertainty of the experimental work. Since 
the magnitude of the variables differs between runs, these values are to 
be considered only as representative. The method of analysis is that pro- 


posed by Kline and McClintock /6/ for single sample data. 


Run Number 3-9-2 2 3-9-10 % 
Item, Units Value + Uncertainty 
m, 1b/min 563.0. 13.6 2.4 534 + 12 22 
Q., ft°/min | 9.03 + 0.22 2.4 | 8.564 0.19 ee 
my 1b 170+ 4 2.4 | 190+ 4 2a 
water 
ty» sec Teel 0.1 0.6 21.4 + 0.1 0.5 
fey: BR 516 + 1 0.2 518 + 1 0.2 
P, psia 18.26 + 0.07 0.4 | 18.26 + 0.07 0.4 
Q,, ft°/min | 0.040 + 0.004 1.0 | 0.325 +0.004 | 0.1 
Oi ét°/min | 0.0562 + 0.006 10.0 | 0.442 + 0.0057 | 1.3 
cae ft°/min | 0.04403 + 0.00439 9.9 | 0.3463 + 0.0047] 1.4 
P, psia 51.9 + 6 10.6 le S 17.5 
ee 0.1121 + 0.0129 10.6 | 0.0368 + 0.0064] 17.5 
Piem? ims He | 30.264 + 0,002 0.007 | 30.264 + 0.002 | 0.007 
en: 76.0 + 0.5 0.65 76.0 + 0.5 0.65 
rere 0.485 + 0.050 10.2 3.89 + 0.77 0.8 
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